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Abstract. Using our new Binary Population and Spectral Synthesis (BPASS) code 
we explore the affect of binary populations on the integrated spectra of galaxies. We 
also explore the interplay of binary populations and a varying maximum stellar mass. 
We compare our synthetic populations to observations of Ha emission from isolated 
clusters and Ha and FUV observations of galaxies. We find that observations tend to 
favour a pure stochastic sampling of the initial mass function although the evidence is 
not significant. We also find that binaries make a stellar population less susceptible to 
the stochastic effects of filling the IMF. Therefore making it more difficult to determine 
if there is a variable maximum stellar mass. 



1. Introduction 

When observations are compared to the results from population synthesis codes perfect 
agreement is rare. There are many uncertain factors that can be tweaked to achieve 
a better agreement. The easiest factor to adjust is the initial mass function (IMF). 
The detail that is never adjusted, beyond altering the metallicity, is the stellar evolution 
models. This is probably because there are only so many sets of available stellar models 
for input into population synthesis and each synthesis code has its favourites. More over 
many stellar evolution groups only provide their best models rather than a collection 
with varying mass-loss, mixing schemes, or binary evolution. 

Creating a full set of stellar models was an intensive process requiring a great deal 
of user and computer time. Today with computer time at a surplus and with recent de- 
velopments leading to more numerically stable evolution codes it is straight forward to 
create many different sets of evolution models with varying input ph ysics. This allowe d 
for large numbers of detailed models of binary stars to be created dEldridge et al.|[2 008). 
Therefore we are now able to vary the uncertain details of stellar evolution within our 
population synthesis code together with other parameters such as the IMF. 

Here we discuss how changing between single star to binary models affects the 
star-formation rate indicators (SFRIs) of Ha and the UV luminosity. We also in- 
vestigate the effect of assuming a varying maximum stellar mass dependent on the 
star cluster mass on these observables. First we outline our numerical method that 
uses the recently developed Binary Population and Spectral Synthesis (BPASS) code 



dEldridge & Stanwa y 2009). Second we discuss how single stars and binary star popu- 



lations differ in their predictions for SFRIs when applied to both star clusters and entire 
galaxies. Finally we summarise our findings. 
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Figure 1 . Schematic diagram showing the inputs and outputs of population syn- 
thesis codes. The boxes with thick borders indicate those unique to BPASS. 



2. Numerical Method 

Recently we have developed a no vel and unique code to prod uce synthetic stellar pop- 
ulations that include binary stars (Eldridge & Stanwayl 120091) . Figure Q] illustrates the 
inputs and outputs that make this code unique. For example the binary models that also 
require other new details such as the initial binary parameter distribution and the effect 
of neutron star and black hole kicks to determine the fate of the binaries after the first 
supernova. We also obtain extra predicted outputs such as the velocity distribution of 
runaway stars. 

While similar codes exi st, for example L eitherer et al.l ( 1999h . iKotulla et all ( 2009 ) 



and IVanbeveren et al.1 (120071) . BPASS has four important features each of which set it 
apart from other codes and in combination make it the state-of-the-art. First, and most 
important, is the inclusion of binary evolution in modelling the stellar populations. The 
general effect of binaries is to cause a population of stars to look bluer at older ages 
than predicted by single-star models. Second is that detailed stellar evolution models 
are used rather than an approximate rapid population synthesis method. Third is the 
use of only theoretical model spectra in our synthesis with as few empirical inputs as 
possible to create a completely syntheti c model to compare with observations. Finally, 



the use of Cloudy (IFerland et all 1 19981) to determine the nebular emission means we 
are modelling not only the stars in detail but also the surrounding gas. 

To provide an example of the importance of binaries in stellar populations we 
show in Figure [2] how a binary stellar population provides a similar spectrum to that of 
a single stellar population with a shallow IMF of -1.35 (where Salpeter IMF is -2.35). 
We see that the differences between a shallow IMF and binary evolution are small, 
on the order of 10%. The difference between them to the normal Salpeter single star 
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Figure 2. Upper panel: the spectra of single star (solid black lines), binary (grey 
line) and shallow IMF (dashed black line) populations. The greatest difference is at 
longer wavelengths. Lower panel: fractional difference between single star spectrum 
and those of a single star population with shallow IMF (solid black line) and binary 
population (solid grey line). 
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spectrum is up to 40%. The main identifying feature of the binary spectrum is that the 
spectrum has more ionizing flux. This leads to the spectral shape being dominated by 
the steps of the nebular continuum emission and the very strong nebular emission lines, 
which may provide a method to distinguish if we are observing a varying IMF or a 
binary dominated population. 



BPASS and its comparison to various observations are discussed in detail in lEldridge et al 



(I2008h and HTdridge & Stanwav (120091) . Some results are available for download at 
http : //www . bpass . org . uk. The code is versatile and adaptable to answer any ques- 
tion that requires knowledge of stellar populations. Here we have altered the code to 
take account of the stochastic nature of the filling of the IMF. We have created multiple 
realisations of clusters and galaxies populating the stellar and cluster IMFs at random 
rather than assuming a smooth IMF or star-formation history. With the code we have 
tested two schemes of filling the IMF. First, a pure stochastic sampling (PSS), where it 
is possible to have a 100M G star in a 100M Q cluster (although it is very unlikely). The 
second scheme is that the maximum mass of a star, M max , in a cluster, of mass M c \, is 
dependent on the total cluster mass. We calculate the maximum star mass from 



log(M max /M Q ) = 2.561og(M cl /M G )(3.82 yi/ + (log(M cl /M G )) yi/ ) OT - 0.38, (1) 



which was derived by iPflamm-Altenburg et al.1 ( 2007 ). In this scheme in a 100M G 



cluster the maximum stellar mass is 9.1M Q . We refer to this as the variable maximum 
mass (VMM) case. In both cases we limit the most massive stars to 120M G . 

We create the stellar population of a synthetic cluster by first setting the maximum 
stellar mass in the cluster from one of the two methods above. I.e. 120M G in the PSS 
case or from the above equation in the VMM case. Then we pick the masses of stars 
in the cluster at random from the IMF. We use a Salpeter slope of -2.35 above 0.5M G 
and -1.3 below this mass to 0.1M Q . We continue to pick stars until the total mass of the 
stars is greater than our target cluster mass. When the last star is added to the cluster we 
examine the difference between the final and target cluster mass, whether the difference 
is smaller with or without the last star. We use the case with the smallest difference. 

The situation for binary populations is slightly more complicated. When we con- 
sider a binary population we assume that each star above 5M Q has a binary companion 
and at random give the star a mass ratio and separation. The mass ratio distribution 
is flat over the range < M%IM\ < 1 and a separation distribution is flat in log a 
over the range of separations from 10 to 10 4 /? Q . We include the mass of the secondary 
companion in the total mass of the stellar population. 

This is the case for a single cluster and can be repeated many times to understand 
the range of possible clusters. We use many synthetic clusters in Section [3J] below. 
When we simulate the stellar population of an entire galaxy we do so by forming multi- 
ple stellar clusters with masses determined from a cluster mass function. The probabil- 
ity of a cluster mass is, P(M C \) oc M~ 2 . We take the minimum cluster mass to be 50M Q 
and a maximum cluster mass of 10 6 M G . A key difference in our method to other such 
simulations is that rather than fixing the star-formation history to some smooth func- 
tion we give each cluster a random age from over 100 Myrs. We then continue to create 
synthetic clusters until the total mass of stars in the galaxy averaged over 100 Myrs 
gives the required star-formation rate. For example to achieve a mean star formation 
rate of lM G /yr requires 10 8 M Q of stars. We find that the random star-formation history 
produces a greater scatter into our predicted observables as seen below. Our reasoning 
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behind this is that nature forms clusters as distinct components and there is no reason 
to expect different clusters to all form at the same time. 

Once we have our synthetic stellar populations of either single clusters with a 
single age or a collection of clusters with a range of ages we can predict various required 
observables. Here we concentrate on two, the Ha and UV fluxes that are both used as 
SFRIs. Ha is provided by stars more massive than 20M Q , while UV flux is produced 
by stars down to 3M Q . Examining these two indicators we have some idea regarding 
how the upper end of the IMF is populated and which of PSS or VMM is in action in 
the Universe. 



3. Results 

Theoretical models are limited in their use without comparison to observations. Com- 
parison both checks the model accuracy and enables us to learn something about the 
physics in action in the Universe. Here we compare our synthetic populations to obser- 
vations of individual star clusters and the emission from nearby galaxies. 

3.1. Individual Clusters 



Recently ICalzetti et al.l ( 2010h have attempted to measure the Ha emission from a num- 



ber of young stellar clusters to determine the mean Ha emission per Solar mass of ma- 
terial in a cluster. A number of clusters have been observed in H a that are selected to 
have ages between 1 to 8Myrs from their broad-band colours. In lCalzetti et all ( 2010h 



the sample of clusters are then binned into two mass ranges. The idea behind this is that 
if the sampling of the IMF is by PSS then a hundred 1000M o clusters should have the 
same stellar population as one 10 5 M o cluster. This means they should have the same 
Ha flux per mass of stars. If there is a VMM then we should expect fewer massive 
star and therefore less Ha per mass of star s. In Figure [3] we compare synthetic clusters 
to the observations of ICalzetti et all (I2010I) . Making sure the models and observations 
have the same mass range and age range. We also include lines that represent the mean 
Ha flux per mass for all the synthetic clusters in the mass range. 

The single star Ha flux decreases more rapidly with age than that for a binary 
population. Thus the single star population does have lower minimum Ha fluxes, by 
about 0.4 dex after 8Myrs. When comparing to the observations there are two pairs of 
points from whether the clusters not detected in Ha are included or not. It is apparent 
that the points in general agree more closely to PSS. Further evidence is supplied by 
the existence of the box in the Figures representing the Velorum cluster that has a total 
mass of 250 to 360M o . It includes two stars in a massive binary which have initial mass 
of 35 and 30M o . For such a cluster VMM predicts a maximum stellar mass of 19M Q . 
While we cannot rule out there is a VMM this one example again points towards PSS 
being more likely. 

3.2. 11HUGS Galaxies 

Another consequence of varying the method of filling the IMF is that different SFRIs 
based on different mass r anges of stars wi ll give discrepant results at different galaxy- 
wide star-formation rates. Lee et al. (2009) estimated the SFR for a sample of galaxies 



from both their Ha and UV fluxes. These results are shown in Figure |4] again compared 
to synthetic stellar populations. In this plot the synthetic galaxies assume constant star- 
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Figure 3. Plot of Ha flux per stellar mass of material in clusters. The points 
are individual realisations of different clusters from two methods of filling the IMF: 
PSS for upper two panels and VMM for the lower two panels. The lines show the 
mean for these two methods, solid black line for PS S and dashed black l ine for 
VMM. The grey crosses represent the observations of ICalzetti et al.l ((2010): light 
grey including clusters not detected in Ha and dark grey not including them. The 
black box represen ts the values derived fro m th e Velorum cluste r estimated from 
Ueffries et al] (120091) . iDe Marco etaf] (120001) and lEMridge (200§- The left panels 
are for clusters including single stars and the right panels are for binary populations. 
The linear features in the points of the upper panels are due to the limited resolution 
of the stellar model's initial masses. 
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Figure 4. The ratio of SFR measured by Ha and UV fluxes v erses the SFR from 
UV flux. The asterisks are the observations from lLee et alj d2009l) while the dots are 
individual realisations of synthetic galaxies. The lines indicate the mean ratios for 
PSS in the solid line and VMM for the dashed line. The left panel is for a single star 
population and the right panel for binary populations. 



formation over a period of lOOMyrs. The models predict that as the SFR decreases 
so to does the ratio of SFRs measured by Ho- and UV. This is because the Ha flux is 
provided primarily by stars of over 20M© while stars down to 3M Q provide UV flux. 
At low SFRs it becomes less likely to get the massive stars required but all less likely 
for them to be formed recentl y enough to contribute to the H a emission. 

It has been suggested by Pflamm-Altenburg et al.l (2009) that the turn down of this 
ratio as evidence for a VMM. However here we show that much of the observed scatter 
can be explained by a stochastic star-formation history. Furthermore this ratio is only an 
upper limit to the expected ratio. There is growing evidence that the covering factor of 
young HII regions is less than unity and therefore many ionizing photons might escape 
galaxies and not produce Ha emission. The UV flux is a direct measure of the flux from 
the stars themselves so is a more reliable measure of the actual mass of stars present. 

In Figure |4]if the single star and binary populations are compared we see that there 
is little difference in the range of scatter. The binary population ratio however does not 
decrease as rapidly as the single star populations. This is because the effect of binaries 
is to increase the number of hot and luminous stars and therefore provide more Ha over 
a longer timescale. Again we sh ould not expect such a rapid turn down as suggested by 
Pflam m-Altenburg et all (l2009h . 



Further work is required in simulating these galaxies but we see here that to deter- 
mine how the IMF is filled cannot be measured from the Ha to UV ratio unless a large 
number of very low star-formation rate galaxies are observed. 



4. Conclusions 

In population synthesis the IMF is an uncertain factor, so to is the evolution of stars. 
Many of the results that are used to infer that the IMF is variable may be explained 
away when binaries are included in population synthesis. 



8 



John J Eldridge 



The Ha flux per Solar mass observed in samples of clusters agrees with PSS of 
the IMF. However more observations are required to rule out a VMM. Individual low- 
mass clusters with one or two massive OB and WR stars are perhaps the most strict 
test to distinguish between PSS and VMM and therefore more such clusters should be 
searched for. 

The ratio of the Ha and UV flux SFRs in galaxies has a significant scatter that can 
be explained by the stochastic nature of the star-formation history, rather than IMF vari- 
ations or a VMM. The ratio of SFRIs for stellar populations including binaries varies 
less than that for populations of single star models only. This is because binaries can ex- 
perience mergers and mass-transfer and produce more massive stars. However finding 
the low-star-formation rate galaxies required to decide between PSS and VMM is diffi- 
cult. Also measuring the SFRIs accurately in such galaxies is difficult. Using a wider 
range of SFRIs and observables beyond those presented here is a more promising av- 
enue to take to discriminate between PSS and VMM. This requires further development 
of BPASS to include more physical processes in greater detail. 
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